A METHOD FOR DETERMINING A RATE OF RAIN 



BACKGROUND OF THE INVENTION 

1 . Field of the Invention 
5 [0001] The subject invention generally relates to a method for 

determining a rate of rain falling on a surface. More specifically, the subject 
invention relates to a method for determining the rate of rain falling on a surface with 
an equation derived from a statistical point process equation. 

2. Description of the Related Art 

10 [0002] Various methods for determining the rain rate are known in the 

art. One example of such a method is disclosed in United States Patent No. 
5,119,002. Specifically, the '002 patent discloses a method that includes the steps of 
detecting a vibration caused by raindrops. The vibration is detected by a sensor which 
produces an output. The method continues by calculating an intermittent period based 

15 on both the sensor output and a constant value. The intermittent period is used to 
actuate a windshield wiper to remove the raindrops. The method of the '002 patent 
does not contain provisions to differentiate vibrations caused by raindrops and 
vibrations caused by non-rain disturbances. Consequentially, the method cannot 
accurately determine a rain rate and will result in improper actuation of the 

20 windshield wiper. 

[0003] Another example is disclosed in United States Patent No. 

5,017,847. The '847 patent discloses a method of detecting the vibrations caused by 
rain, converting the vibrations into electrical signals, then using the electrical signals 
to directly actuate the windshield wiper. 
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[0004] Due to the deficiencies in the methods of the prior art, 

including those described above, there remains an opportunity to introduce a method 
that is capable of determining a rain rate based on the principles of a point process. 
The point process provides a means to determine a very accurate rain rate, since rain 
5 exhibits characteristics of a point process. Conversely, non-rain disturbances are 
inherently removed from determination of the rain rate, since non-rain disturbances do 
not exhibit the characteristics of a point process. 



SUMMARY OF THE INVENTION AND ADVANTAGES 
10 [0005] The invention provides a method for determining a rate of rain 

falling on a surface. As rain falls and strikes the surface, the surface vibrates. The 
method includes the step of sensing the vibrations of the surface. A vibration signal is 
generated, the vibration signal being proportional to the vibrations of the surface. The 
vibration signal includes peaks. The method also includes the step of determining the 
15 peaks of the vibration signal. Time intervals between the peaks are then determined. 
Using the time intervals, a number ni of time intervals that occur between a first time 
and a second time are counted. The number ni of time intervals falls within a first 
range. A number n2 of time intervals that occur between the second time and a third 
time are also counted. The number n 2 of time intervals falls within a second range. 
20 The rain rate X is then determined using an equation that is derived from a point 
process equation and utilizing ni and n 2 . 

[0006] The natural vibrations caused by rain exhibit the properties of a 

point process. However, vibrations caused by non-rain disturbances do not exhibit the 
properties of a point process. An accurate rain rate is established by the subject 
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invention due to the use of the equation derived from the point process equation. The 
point process equation inherently differentiates between vibrations caused by rain and 
vibrations caused by non-rain disturbances, such as audio waves, wind, and a running 
engine. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0007] Other advantages of the present invention will be readily 

appreciated as the same becomes better understood by reference to the following 
detailed description when considered in connection with the accompanying drawings 
wherein: 

[0008] Figure 1 is block diagram illustrating a method for determining 

a rate of rain falling on a surface according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
15 [0009] Referring to Figure 1, wherein like numerals indicate like or 

corresponding steps, a method for determining a rate of rain falling on a surface is 
generally disclosed at 10. As rain falls and strikes the surface, the rain causes the 
surface to vibrate. Preferably, the surface is a glazing of a vehicle, but could be any 
surface capable of vibrating when struck by rain, such as a hood or roof of a vehicle. 
20 Those skilled in the art will appreciate that the term "glazing" refers to a windshield, 
rear window, side window, etc. of a vehicle. Other disturbances from non-rain 
sources, such as wind, audio waves, a running engine, etc., may also cause the surface 
to vibrate. The method described herein determines a rate of rain falling on the 
surface while inherently removing disturbances caused by non-rain sources. 
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[0010] Referring to Figure 1, one embodiment 10 of the method of the 

present invention begins with the step 12 of sensing the vibrations of the surface. 
Preferably, a piezoelectric vibration sensor is connected to the surface for sensing the 
vibrations of the surface. However, those skilled in the art will realize that other 
5 sensors or apparatus for sensing the vibrations could also be utilized, such as 
potentiometric or variable reluctance-type accelerometers. 

[0011] The method continues with the step 14 of generating a vibration 

signal proportional to the vibrations of the surface. The vibration signal includes 
peaks representing vibrations that occur when a raindrop strikes the surface. The 
10 piezoelectric sensor of the preferred embodiment generates an analog electrical signal 
as the vibration signal. Other non-rain disturbances may also cause peaks in the 
vibration signal. 

[0012] The method also includes the step 16 of determining the peaks 

of the vibration signal and the step 18 of determining time intervals between the 
15 peaks. It is preferred that these steps are performed by a microprocessor that is 
operatively connected to the piezoelectric vibration sensor. However, those skilled in 
the art will appreciate that other devices, including but not limited to microcontrollers, 
application specific integrated circuits, and analog circuits, are capable of performing 
these steps. 

20 [0013] The method continues with the step 20 of counting a number ni 

of time intervals that occur between a first time and a second time and that fall in a 
first range of the time intervals. In the next step 22, the method counts a number n 2 of 
time intervals that occur between the second time and a third time and that fall in a 
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second range of the time intervals. Those skilled in the art will appreciate that the 
first and second ranges are also known as "statistical bins", or simply "bins". 

[0014] The rain rate X is determined with an equation that is derived 

from a statistical point process equation in step 24. An exponential density function 

of a first order point process is represented by the equation f(t) = Xe~ M , for t > 0, 
where X represents the rain rate, t represents time, and e represents the exponential 
function. Those skilled in the art understand that the exponential function e is an 
irrational number whose numerical value is approximately 2.71828. The vibrations 
caused by rain exhibit the characteristics of a point process, also known as a Poisson 
process. However, other vibrations caused by non-rain disturbances, such as audio 
waves, wind, and the running engine, do not exhibit the characteristics of the point 
process. Therefore, the vibrations caused by the non-rain disturbances do not 
significantly effect the determination of the rain rate X, according to the present 
invention. 

[0015] In a preferred embodiment, an analog-to-digital converter 

(ADC) is operatively connected to the vibration sensor. The microprocessor is 
operatively connected to the ADC. The method further includes the step of 
generating a series of digital values that are representative of the vibration signal. 
Each digital value has a magnitude component and a time component. Preferably, the 
ADC generates the magnitude component and the microprocessor generates the time 
component. However, it should be appreciated that different hardware configurations 
can be used to generate the magnitude and time component of each digital value. 

[0016] The method of the preferred embodiment also includes the step 

of modifying the magnitude components of the digital values to enhance the peaks of 
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the vibration signal. Enhancing the peaks of the vibration signal allows for easier 
determination of the peaks by the microprocessor. The step of modifying the 
magnitude components of the digital values to enhance the peaks of the vibration 
signal is preferably broken down into three separate steps. First, the step of detecting 
a first magnitude v(n-l), a second magnitude v(n), and a third magnitude v(n+l) is 
performed. These three magnitudes are consecutive digital values. Next, an 
enhanced magnitude v e (n) is computed with an equation v e (n) = v 2 (n) - Av(n- 
l)-v(n+l), wherein A is a scaling number selected to maximize the enhancement of 
peaks. Finally, the step of successively determining enhanced magnitudes v e (n) to 
provide an enhanced series of digital values is executed. Each enhanced digital value 
has an enhanced magnitude component and the time component. 

[0017] In order to determine the peaks of the vibration signal, the 

preferred embodiment uses a multi-step process. Preferably, this process is performed 
by the microprocessor. The process begins with the step of sorting the enhanced 
series of digital values by the magnitude component of each digital value. The 
process continues with the step of discarding all digital values with magnitude 
components that are not greater than a certain threshold from the enhanced series of 
digital values. The certain threshold is determined to separate signal "noise" from 
significant disturbances to the signal. The certain threshold is based on various 
factors, including, but not limited to, sensitivity of the vibration sensor, amplification 
of the vibration signal, and resolution of the ADC. An ancillary step in the process is 
determining the certain threshold based on a size of a droplet of rain. For example, in 
a rainstorm with small and light droplets, the certain threshold would be different 
from the certain threshold in a rainstorm with large and heavy droplets. 
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[0018] The process continues with determining the largest enhanced 

magnitude in the enhanced series of digital values to identify a valid peak. Then, a 
next largest enhanced magnitude component in the enhanced series of digital values is 
determined to identify another valid peak. The next largest enhanced magnitude is a 
5 set time apart from the largest enhanced magnitude component. The set time, 
typically in the order of milliseconds, is based on a number of factors (resolution of 
the ADC, etc.). The process concludes with the step of successively determining the 
next largest enhanced magnitude component in the enhanced series of digital values 
that is a set time apart from any of the valid peaks to identify additional valid peaks. 

10 Once a valid peak is determined, the process ensures that another peak is not found to 
be "valid" within the set time on either side of the valid peak. This assures that up 
and down swings related to a peak are not identified as valid peaks. 

[0019] In another embodiment, a different technique is used to 

determine the peaks of the vibration signal. This technique begins with the step of 

15 analyzing each digital value as the digital value is generated from the digital signal. 
This is done to determine if the magnitude component of the digital value exceeds a 
certain threshold. This technique can also be practiced by substituting the enhanced 
digital value for the digital value. The technique continues with the step identifying a 
valid peak when the magnitude component of the digital value first exceeds the 

20 certain threshold. The technique concludes by starting a timer to run for a certain time 
when the magnitude component of the digital value exceeds the certain threshold. 
The technique does not allow the identification of other valid peaks while the timer is 
running. The length of the timer is longer than a refractory time caused by a droplet 
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of rain striking the surface. Also, the certain threshold may be varied based on a size 
of droplets of rain. 

[0020] In order to easily determine the rain rate X using the 

microprocessor, three different alternative embodiments have been developed. 
However, it will be realized, by those skilled in the art, that other alternative 
embodiments based on the point process could also be developed. 

[0021] A first alternative embodiment involves solving for the rain rate 

X by taking a mathematical derivative of the exponential density function of a first 
order point process. 

[0022] The derivative of f(t) is f{t) = -AV* , which can be rewritten 

as f y (t) = -Z- f(t). Solving for the rain rate X yields A = - • 

[0023] The equation for the rain rate listed above used functions that 

are continuous over time. The microprocessor of the preferred embodiment utilizes 
discrete values to find the rain rate X. Therefore, the equation is 

rewritten X = — , where m represents the ranges, or statistical bins, of time 

intervals that are used to discretely approximate f(t) and f (t). 
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[0024] Two of the statistical bins can then be used to determine an 

intermediate rain rate 7^ nt . By solving for f (m) and f(m) 5 the equation for determining 
the intermediate rain rate Xj n t is determined, where ni and n 2 are the statistical bins and 
w is the width of the statistical bin having a unit of time. 

f{m) = r hz2h_ 
w 

f( m ) = 2 2 

2 n 2 - n x 

Am ; — 

w n 2 + n x 

Using several of the statistical bins, multiple solutions for Xi nt are determined. 
Averaging the multiple solutions for hm together yields the rain rate X. 



[0025] The first alternative embodiment for solving for the rain rate X 

can be summarized as follows. First is the step 46 of maintaining the first range and 
the second range of time intervals equal in a time span w. All statistical bins must 
have the same size with respect to time. Next, the step 48 of determining an 

2 J2 

intermediate rain rate ^ nt with an equation X mt = L is performed. The 

w n 2 + n x 

method of the first alternative embodiment continues with the step 50 of successively 
determining the intermediate rain rate Xi nt . For example, in a first iteration a first and 
a second bin are used to determine the intermediate rain rate A^ nt . In a second 
interation, the second bin and a third bin are used to determine another intermediate 
rain rate A^ nt , and so on. Finally, the intermediate rain rates A^t are averaged together 
to determine the rain rate X in the concluding step 52. 
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[0026] A second alternative embodiment involves solving for the rain 

rate X by using an approximation with only two statistical bins. The first range of the 
time intervals is defined as all time intervals less than or equal to the second time. 
The second range of the time intervals is further defined as all time intervals greater 
5 than the second time. 

[0027] To show how the equation of the second alternative 

embodiment is determined, first, suppose the exponential density function of a first 
order point process is rewritten in discrete form / (m) = Ae'^ . 



The first bin is then estimated 



/(i) = ^ 



-X 



m=2 m=\ 

10 The second bin can be estimated as 
Now solving for the rain rate 



x = /(l) = tor* _ n x 



e n 2 



m-2 



[0028] The second alternative embodiment includes the step 54 of 

determining the rain rate X with an equation X = ri\ I n2. 
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[0029] 



A third alternative embodiment also involves solving for the 



rain rate X using only two statistical bins. However, an equation of the third 
alternative embodiment will provide a more accurate determination of the rain rate X 
than the equation of the second alternative embodiment. The first range of the time 
intervals is defined as all time intervals less than or equal to the second time. The 
second range of the time intervals is further defined as all time intervals greater than 
the second time. 

[0030] The equation of the third alternative embodiment is determined 

by starting with the equation of the exponential density function of a first order point 




process 



[0031] 



and understanding that 



00 



\f{t)dt = 1 



0 



[0032] 



Then, 



T 



T 




0 



0 



where T represents the second time and "a" is a constant. Continuing, 




T 



0 



where t is substituted by 9 +T. 
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[0033] 



Now solving for the rain rate X, 



l-e 



e 



^ = e XT -\ 




In -^ + 1 = 



J 



\ 



when AT is small. 



[0034] 



The third alternative embodiment includes the step 55 of 



determining the rain rate X with an equation X = ni / (T * n 2 ), where T represents the 
second time. 



further includes the step 56 of producing a rain signal in response to the determined 
rain rate X. Typically, the rain signal is a digital signal that is transmitted to a wiper 
computer. Those skilled in the art will appreciate that the rain signal could also be an 
analog signal. 



wiper computer then uses the rain signal to control a frequency of a windshield wiper 
based on the rain rate X. Furthermore, an additional step of modifying the frequency 
of the windshield wiper is performed. Modification of the frequency is performed 
utilizing a hysteresis process. The hysteresis process uses the current frequency of the 
windshield wiper and the rain rate to modify the frequency, resulting in a smooth 
transition of the frequency. 



[0035] 



After the rain rate X is determined, the preferred embodiment 



[0036] 



Once the rain signal is received by the wiper computer, the 
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[0037] Obviously, many modifications and variations of the present 

invention are possible in light of the above teachings. The invention may be practiced 
otherwise than as specifically described within the scope of the appended claims. 
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